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SBCTIOR 1.0 
INTRODUCTIOil 


Laser tracking of the LAGBOS satellite la presently 
being perforaed with noise levels at the few centlaeter 
level and with systeaatlc errors (biases, refraction, 
tlalng, etc*) considered to be at or reducible to the sub<> 
centlaeter level. To take full advantage of such 
accuracies, the aeasureaent and epheaerls aodellng In the 
data reduction program must have comparable accuracies. The 
aost widely used computer prograa for laser data processing 
Is the 6SFC Geodyn prograa* In Its original foraulatlon, 
this prograa Integrates satellite equations of aotlon based 
on Newtonian mechanics, and allows the estimation of orbital 
and geodetic paraaeters so that a best fit Is obtained. In a 
weighted least squares sense, to an Input data set 
consisting of various data types. This report describes the 
modifications to Geodyn to substitute the Einstein 
gravitational theory for Newtonian gravitation. This results 
in aodlf tcatlons to both the satellite equations of aotlon 
and to the aodellng of satellite tracking aeasureaents . 


Since the Newtonian theory Is a very close 
approxiaat ion to the Einstein theory (or general theory of 
relativity), observable deviations from Newtonian theory are 
small, and the applicable equations of aotlon or observation 
equations can be formulated so that they differ only by 
small terms which we will denote as relativistic 
perturbations. In the implementation of the Einstein theory 
In an orbital data reduction prograa, one still deals with 
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station coordinates and satellite coordinates which resewhle 
the noraal Cartesian coordinates* one has a tlae systea 
defined for the coordinate systea* and one still aeasures 
round trip transit tlaes froa ground stations to 
satellites. However* the satellite equations of aotlon are 
altered though the addition of coapllcated perturbation 
terns* the clocks at tracking stations do not aeasure the 
saae tlae used In the equations of aotlon* and one does not 
obtain a range aeasureaent by slaply aultlplylng the transit 
tine by a constant. And one aust continually be aware that 
paraaeters considered constant In the Newtonian theory aay 
no longer be considered so In the relativistic theory. 

In general* no atteapt will be aade In this report to 
derive the perturbation equations which should be 
Inpleaented t ' transform Geodyn (or comparable computer 
program developed for estimating orbits for earth 
satellites) into a relativistic program. Mout of the 
desired relations are well documented In the literature and* 
in some cases* have a very complex derivation. Likewise* no 
extensive discussion will be given of general relativity 
theory Itself, for which textbooks aay be consulted [e.g.* 
Misner, Thorne and Wheeler, 1973]. The basic theory 
deviates from classical mechanics in that the "Interval" 
between points is the Integral /ds* where 


ds 


u V 

g dx^ dx 


( 1 ) 
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is the metric tensor 
W, V - 1 , 2, 3, 4 

and summation of |i,v from 1 to 4 in (1) is implied* The 
coordinates Include, of course, three spatial coordinates 
and the time coordinate. The path followed by any 
unconstrained body is such that the Integral /ds is a 
maximum. In Newtonian theory, a coordinate system can be 
chosen in which g ■ *1* 5 and one has Cartesian 

coordinates* In the Einstein theory, the g^^ depend upon 
the distribution of mass bodies, but a coordinate system can 
be chosen for the motion of planetary bodies such 

chat g^^ is almost the Minkowski metric (diagonal 

-1 ,-l ,-l ,+l ) . In fact, the only available solutions for the 
motion of planetary masses are in such a system* 

In the implementation of relativistic effects, two 
Important points must be emphasized: 

1* The coordinate system for integration and 
measurement computation becomes a solar system 
barycenter system. The equations of motion for 
the satellite are derived in this system, and it 
has already been implicitly assumed in the 
gravitational third body perturbations in orbit 

programs such as Geodyn* In practice, Geodyn 
normally Integrates the satellite acceleration 

relative to the center of the earth in order to 
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nalntain highar praeiaton. Cara thus naada to ba 
axarclaad In tha obaarvatlon eoaputatlon to 
inaura that, aquivalantly , coordlnataa rafarancad 
to tha aolar ayatan barycentar, or diffarancsa of 
auch coordlnataa, ara baing uaad* 

2. A nuabar of ehangaa are required to laplaaent 
Elnateln gravitation Into a prograa auch aa 
Geodyn. The net affecta of all ehangaa on 

obaervablaa (auch aa baaellnea between atatlona) 
are expec.ed to be aaall (e.g., below tha 
decimeter level). This doea not neceasarlly mean 
that the effects of Individual changes are 
small. The relativistic theory Is different from 
Newtonian theory, but must be accepted only as a 
complete theory. 

This report Includes consideration of both of these 
points. The effects of using an explicit earth-centered 
coordinate system will be examined with regard to the 
computation of the laser ranging observable. And simulation 
results will be presented to demonstrate the effects of the 
Individual modifications required to Implement the Einstein 
gravitation theory In an orbital data reduction program 
designed for earth orbiting satellites. 
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SECTION 2.0 

RELATIVISTIC SOLUTIONS FOR n-BODT MOTION 


Although quite elegant in teneor forauletlon» the 
Einstein equations for determining the ten metric tensor 
components are quite difficult to solve for a general set of 
physical bodies* For sate Lite motion, the solution needed 
is one which can account for the motion of the satellite in 
the gravitational field of the earth, sun, and moon. The 
influence of the other solar system bodies should be 
considered, but relativistic effects from even Jupiter and 
Saturn would not be expected to be significant in comparison 
with those of the sun and earth* 

The solution of the Einstein field equations for a 
single point mass was first obtained by Schwarsachild [1916] 
within a year of the publication of the general relativity 
theory* Shortly thereafter, approximate solutions for the 
field of n mass points [Droste, 1916] and for the motion of 
n heavy bodies [de Sitter, 1916; Eddington and Clark, 1938] 
were also obtained* The approximate n**body solution should 
be perfectly adequate for the application to earth satellite 
motion* The most appealing approach to the motion problem 
has been that taken in the so-called Einsteln-lnf eld- 
Hoffmann (ExH' approximation technique, since the equations 
of motion wer<i obtained directly from the field equations* 
These equations of motion give the accelerations for a 
spherically symmetric body moving in the field of n-1 other 
spherically symmetric bodies, in a non-rotating coordinate 
system centered at the center-of-mass of the n bodies* To 
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b« coniittcnt with th« tolution dwrlwwtlont thtsw 
Aecwlwrwtioni auit bw wffwetiwwly eowputwd la tha aolar 
•yat«B barycantar ayataa for both tha aatalllta aod for tha 
aarch. Station eoordlnataa suit alao ba In tha aa»a ayata«. 

Tha coordlnata ayatan ganarally adoptad for tha n- 
body problaa haa baan a non-rotating ayataai whoaa o*'lgln la 
at tha cantar of aaaa of tha n-bodlaa* Tha eoordlnataa In 
thla ayatan ara naarly Cartaalan and tha aquatlona of motion 
for tha bodlaa ara almllar to thoaa for Nawtonlan n-body 
motion but with amall additional ta&.'^a. In tha Gaodyn 
implamantatlon, thaaa additional tarma ara traatad aa tha 
ralatlvlatlc parturbatlona, with tha Intagration otharwlaa 
procaadlng aa for Intagration of tha Nawtonlan aquatlona of 
motion. Abaoluta tlma in tha Intagration la raplacad by 
coordlnata tlma. 

Tha accalaratlon of one of tha n-bodlaa In thla 
canter of aaaa ayatam may ba obtained by carrying out tha 
coordinate timn differentiation in tha EIH aolution (Infald 
and PlabaiTakl, p. 112]. In the notation of Moyer [1971, Eq. 
35], the raault after delating the Newtonian term la: 
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opraoRQUMirr 


if. 


J7i r:, i#i lA c* kJj 'jk 


'ij 


* at* ‘**‘*‘’ ■ '*1 * t* ^ 


. _z_ 2 : ijji 

2c* jji 'ij 


( 2 ) 


wh«r« 


x»y»* 

• • k 

•• •• •« 
«»y.* 


gr&vitittlon«l constant for body J* 


rsctsngulsr conponsnts of position, 
velocity, snd sccslerstion relstlvs to s 
non-rotsting frame of reference centered at 
the barycenter of the system of n bodies. 


r^j - l(xi-Xj)^+(yi-yj)^+(*l-8j)^J 


1/2 


.2 ^ .2 ^ .2 

*1 + yi + 


All dots denote differentiation with respect to coordinate 
time • 
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This accslaratlon is ths rslstlvlstlc perturbation of 
a body, and should bs valid for both the aceslsration of ths 
LACBOS satallits in ths solar systsa barycsntsr systsa, and 
ths aceslsration of ths sarth in ths saas systsa. Ths 
rslativisbically coaputsd aceslsration of ths diffsrsnes 
bstwssn ths satsllits and sarth coordinatss is thus 
obtainabls by differencing the satsllits acceleration froa 
(2) and ths sarth aeeslsraticn froa (2), or, 




(3) 


Ws call this diffsrsnes ths acceleration of ths satsllits 
relative to ths earth. Ths satsllits-sarth coordinate 
differences aay, using these perturbations, bs integrated as 
an alternative to separate intsgrstion of satsllits 
coordinates and sarth coordinates. A considerable increase 
in precision is thereby achieved. In practice for the 
LAGEOS satellite, only tho sun and sarth constitute 
significant perturbations in Eqn. (2) or (3). 

Associated with the n-boly equations of notion Is a 
netrlc tensor and line elenent. Again in the notation of 
Moyer [1971, Eqns. 22-30], the n-body l«ne eleaent can be 
written 
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1*1 'ij ^ 


(dO* 



PH 

( j#l ’'ij c* 

Ijj 


J#1 'ij 


c^ J*i ij k*j 


y;„ [!!£Vi!i * iL . iiiciiiiii 

c* f^i J I '11 '11 r\ 


'ii 




(4) 


This line element will be needed In the v^erivetlon of the 
relativistic range calculations. It could also be used in 
developing the equations of smtlon of the i'th particle by 
finding the trajectory which extremizes /ds. The result is 
Eqn. (2) (plus the Newtonian effects)* 
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SECTION i.O 

RANGING MEASUREMENT CALCULATIONS 


The baaic ■aasureaent aade by a aatalllta ranging 
station la the time that elapses on the station clock 
between the transmission and reception of a signal* For the 
laser tracking case* the signal Is simply reflected, so the 
delay at the spacecraft* Is relatively easily calculated, at 
least for a spherical satellite suc'i as LAGEOS. The problem 
Is then to relate the (corrected) time interval measurement 
to the coordinates of the tracking station (at transmission 
and reception) and the satellite* 

This problem has been addressed a number of times for 
the most significant perturbing body, with the result [e*g*, 
Moyer, 1971, p. 171 


^11 

'j - - c + ^ 

•* c 


^1 ^Ij 


'l ■ "iJ 


(5) 


with the signal traveling from point 1 at coordinate time tj. 
to point J at coordinate time tj. r^ and rj are the 
magnitudes of the position vectors from points 1 and J to 
the sun and Is the gravitational constant of the sun. 

d 


*Delay here means the correction necessary to the measurement to 
make It effectively a measurement to the spacecraft center of 
mass • 
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This r«tult can alto be obtained ualng the line 
element from Eqn. (4). If we retain only firet order 
perturbation terms (for an earth orbiting satellite), drop 
the subscript 1, and consider that the path described Is a 
null geodesic, we have 



where 


r 


J 


4 


x-x 


j 






and J ranges over all perturbing bodies* If all body motion 
la conaldered to be slow, or to be negligible compared to 
the total transit time, then we can set 

X “ Xj + o (X2 ~ Xj) 

y - Yj + ® 

z ■ Z2 + ^ ^^2 ~ 
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where 


(Xj, Y^. Zj), (X^, Y^, Zj) 


are the coordinates of the transaiitter and receiver, 
respectively. The paraaeter a has the range 0-1. 

Substituting Into Eqn. (6), transposing the first term to 
the left, and taking the square root of both sides, we 
obtain, approxlaately 



Integrating both sides, 

^ *’12 do 

c At - r., + — 2 — " 7“ 
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QIMUIV 


To perfora the last Integration* we write 


r^ - /(x-Xj)^ + (y-y^)^ + (a-a^)^ 


- / 


[Xj+a(Xj-X^ )- l^+[ Yj+odj-Yj )-y^ J^+CZ^+aCZ^-Z^ ^"*J 


j j_ 

2 “ ^2 • “ ^12 


( 10 ) 


Substituting Into (9) and performing the Integration* 


c 



*^2 ■*■ ^12 *12 


*J1^^12 


r 4- g 

*^1 *12 




The argument of the logarithm can be shown to be equivalent 
to the argument of the logarithm in Eqn. (3)* so 
that c At can be written 


c At 


^12 ■*■ 




In 


Ti + r2 + rj2 


'i “ ’^IZ 


(11) 
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The total relativistic perturbation is thus equal 
of perturbations of the Individual effects 

gravitational bodies. The largest effect for 

orbiting satellite such as LA6E0S is the sun» but 
is also significant at the centlneter level. 


to the sun 
of the 
an earth 
the earth 
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SECTION 4.0 

STATION COORDINATE CORRECTIONS 


The contraction of aovlng rods la veil known fro« 
special relativity theory. In general relativity, a alallar 
phenomenon occurs except that the term "rigid** must be 
defined. The subject of solar system barycenter coordinates 
for tracking stations fixed on a rigid earth* has been 
considered by Nlsner [1982]. using the concept that the 
proper length /ds from the earth center world line to the 
tracking station world line remains constant, with the 
distance measured In a spacetime direction orthogonal to the 
earth's world line. This concept provides a scalar 
relationship between the solar system barycenter coordinates 
(which must be used in the data reduction program) and the 
station coordinates moving with the earth. On intuitive 
grounds, along with this scalar condition. Mlsner deduced 
the following relation between the two coordinate systems: 


- X- - [U„ + (U_ 1/2 V_^) 1 X_/c^ 

T TEE E avg T 


- I • *-f> 


( 12 ) 


*Dynamic motion of stations Is assumed to be otherwise 

accounted for. 


15 


where 


X«, are the station coordinates in a coordinate 

X 

systea aovlng with the earth but non-rotating* 

are the station coordinates in the solar systea 
frame. 

▼ is the earth velocity in the solar systea frame. 

£ 

U£ is the gravitational potential due to the sun at 
the earth. 

The net effect is a nearly constant scaling of the station 
coordinates by 2 1/2 parts in 10^ and a dally variation 
which can be as large as ^ 3 cm. 

It should be emphasized that, while Eqn. (12) may be 
intuitively appealing and may be used with a reasonable 
degree of confidence, it has not yet been rigorously derived 
or verified. 
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SECTION 5.0 

COORDINATE-ATOMIC TIME TRANSFORMATIONS 


By definition, an interval of "proper” tiae aeaaured 
by an atoaic clock la proportional to the interval da along 
its world line, 


dT ■ 



(13) 


The clock time interval is thus not proportional to the 
interval of coordinate tine which is used for trajectory 
integration. However, the scale of the clock is, by 
convention, deliberately chosen so that atomic time and 
coordinate time are, on the average, colinear. The n-body 
line element given by Eqn. (4) can be used for ds in Eqn. 
(13), and the deviations between t and t can then be 
calculated on the basis of the motion of the clock and the 
ephemerides of the gravitational bodies. 

The coordinate time - Atomic time differences (t-x) 
have been analyzed by various investigators, in particular 
by Moyer in two recent papers [Moyer: 1981a, 1981b]. 

Moyer's final result, in a form suitable for Implementation 
in a computer program such as Geodyn in which the planetary 
ephemerides are conveniently available, is: 
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t - 1 . «A + (i® . r|) + ij <ij . ^ (ij • *;> + 


c*(li„ + u.) 




SA 


(i 


S 


2, ^ V "SA 


8 \ « 1 ^ • C S \ 

'SA^ * 2 

c 


(14) 


The ■ubecrlpts and superacrlpte In this equation refer to: 

A ■ location of atonic clock on Earth which 
reads International Atonic Tine t 


E 


Earth 


Earth-Moon barycenter 


M 


Moon 


Sun 


Solar aysten barycenter 


Jupiter 


SA 


Saturn 


The quantity dT^ Is a constant offset (32.184 seconds), the 

p's are the nasses of the subscript body, and 

t' and r^ (e.g.) denote the position and velocity of the 
B B 
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•«rth-aoon baryc«nt«r (aubtcrlpt) with rctpact to tho ion 
(tupcrocrlpt ) • 

Tho largoat amplituda non-conatant tor* on tho right 

hand aido of (14) la the torn 2(f| • r|)/e^f with an 

anplituda of 1.658 aaec and a 1 year period. Teraa due to 

Jupiter and Saturn have aaxiaua awplltudea on the order of 

20 ysec and perloda allghtly greater than 1 year. The only 

• C E 2 

abort period term la <r, • r.)/e , which haa a maximum 

E A 

amplitude of ~ 2 paec and a period of 1 day. 

It should be noted that the expreaalon for t-t above 
is approximate, and with the retained terms chosen on the 
basis of their Influence on the NASA Deep Space Network 
observables (range and range rate). However, this criterion 
resulted in the retention of terms with amplitudes greater 
than 3.7 usee, 0.11 Msec, and 1.3 Msec for periods of a day, 
a month, and a year, respectively [Moyer, 1981a]. Eqn. (12) 
should thus be more than adequate for the processing of 
earth-orbiting satellite data. 
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SECTION 6.0 

MEASUREMENT COMPUTATION IN 
BARTCENTER COORDINATES 

At hat bttn indicated, the utilisation of the EIH 
tolution for the tatellite equetiont of aotion hat iaplied 
the adoption of a tolar tpttea barycenter coordinate 
tyttem. In terae of the light tiae tolution given in 
Section 3.0, thit aeent that the difference between pulte 
trentaittion end return tiae at a tracking ttation it 
expretted at: 



(15) 

where 

tf it the later firing tiae 

tjj it the later receive ("bounce") tiae at the 
tatellite 

tp It the receive tiae at the later tite 
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r| it tht rang* ttom tht flrlog point to tbt tun at 
tint tf 

r| it tht racga froa tht tatallitt to tht tun at 
tint t^ 

r| it tht ranga froa tht rtcaiving tita to tha tun 

at tint tp 

^ it tht ranga froa tha trantaitttr point at tiaa 

t£ to tha tatallitt at tiaa t^^ 

it tha ranga froa tha tatallitt at tha tiaa t|^ to 
tha rtcaiving aita at tiaa 

For aia;*licity, ralativiatic affactt of bodiat othar than 
tha tun hava baan naglactad in Eqn. (15). To ralata tf^ + 
r|jp to tha ranga calculation noraally parforaad in coaputar 
prograat operating in earth"‘canterad coordinatati wa write 
out thaaa rangaa explicitly and expand all eoordinataa about 
their valuta at tha tatallitt tiaa t^^, Firtt, wa denote 

Xl» Zi at the tatallitt eoordinataa 

Xt» Yg, Zg at tha tracking ttation eoordinataa 

Xgi Ygt Z£ at tha eoordinataa of tha canter 

of tha earth 

All theta eoordinataa are in tha tolar ayttem barycantar 
tyttea. Tha range tun tf^ + r^^^ it than defined by 
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'fb * 'br 


. /r. 


*l<'b’‘*.<*f‘'* * I’'L<*b>-*,<‘f>l* + l*L<‘b>-*,<'f>l’ 


♦ Aj 




( 16 ) 


Next consider the coordinate difference 

• VS>-*b“b> - l*,«'b>-*E“b” 

- (t, - t^) - . . . 

‘ W'b> - *B."b> • */*b’ ‘'f - S> - • • • 

* - *.‘'b’ “f * ‘b> * • • • 
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Haking uta of (17) In (16) » tha lattar baconaa 
'fb ♦ 'bt ■ /* - J « • + <*.>* 

+ /r - 2 t . i,(t^-t^) (i^)2 

(18) 


whare 


• ■ ‘*.l"b>- ’'.L<S>* ^L"b» 


*. ■ *.<‘b> 


Expanding Che square roots In (18) » we obtain 


R • X 


""fb ’'br ^ - ‘k + (t. “ K)1 


f b '"r "b 


1 ^*s> 2 2 

~ 1 


( 19 ) 
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To first order we can set 







2 (R/c)^ 


Eqn. (19) then reduces to 



( 20 ) 


Substituting this result Into Eqn. (15) and making minor 
approximations In the logarithmic term» we obtain 



I (21) 

from station and satellite, 
evaluated at the satellite 

time t^« 


2 GM, 


in 


r® + r? + . 
s L 

S ^ S 
r + r, - I 
s L 


s s 

where r and r. are the rang 

8 L 

respectively, to the sun, a 
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Except that it It beeed on coordinate rather than 
atoaic tiaaa, the left hand aide of (21) ia the noraal 
aaasurad range and the first tera on the right hand aide ia 
the noraal calculated range. The second tera on the right 
hand side ia a correction required hacauea of the uaa of 

earth centered coordinates and the coaputation of a range 

baaed on an average set of station coordinates. The last 
tern on the right hand side ia a true relativistic 

correction. 

It is iaportant to recognise that the predoainant 

part of in the second term arises froa the earth's 

barycentric velocity of 30 ka/sec» and is only slightly 
modified by the earth's rotational velocity of ~ 400 a/sec. 


SECTION 7.0 
SIMULATION RESULTS 


To assess the aagnlttide of the relativistic 
perturbations In the reduction of LAGEOS laser tracking 
data, simulations were performed using a modified version of 
the ORAN error analysis program. With the exception that 
the range measurement corrections from Sections 3 and 6 were 
lumped together, all the perturbations were carried 
separately. The total effects on the estimated parameters 
are then the algebraic sun of the Individual effects. 

The simulations considered the following data set: 

Laser Stations 

- Goddard Space Flight Center (STALAS) 

- Arequlpa, Peru (ARELAS) 

Yaragadee, Australia (YARLAS) 

Owens Valley, California (OWENSV) 

Wetzell, West Germany (WETZEL) 

- Haleakala, Maul, Hawaii (HOLLAS) 

Station Visibility 

- 50Z of passes tracked 

Elevation Cutoff - 20^ 

Arc Lengths: 2 weeks 

Number of Arcs: 6, spaced two months apart In 1980 
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Two types of siaulatlons were performed. First, 
noting that data reductions with the relativistic 
perturbations will require a different scaling fron GMg, all 
6 arcs were processed with the conaon estiaatlon of 
coordinates for all stations (with one station longitude 
held fixed), the coaaon estiaatlon of GMg, plus the 
estiaatlon of the 6 orbital elements for each arc. The 
results of this simulation are shown In Table 1. For the 
station positions, the dominant effects are those of the 
station coordinate modifications and measurement modeling on 
the station heights. These effects, however, almost cancel, 
leaving a net height effect on the order of a centimeter. 
For the estimated GM^^ the force model and measurement 
modeling produce the only significant effects. There is 
about 751 cancellation here, leaving a net change In GMg of 
-.00824 km^/sec^. Since the sign convention in Table 1 Is 
based on applying relativistic corrections versus not 
applying them, reduction of the simulated data set with the 
relativistic models will result in the estimation of a 
smaller value of GM^ than would be estimated using pure 
Newtonian models. 

We next consider that the 6 arcs are used 
Independently to estimate the same set of station 
coordinates, with a fixed value of GIl£. For comparison with 
a Newtonian solution, we assume that the relativistic 
solution uses a value of GM£ which is .00824 km^/sec^ 
smaller than that used in the Newtonian solution. Table 2 
shows the results for one station for one of the 6 arcs. 
There are effects from all perturbations, except for timing. 
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station Force Model Station A1*ET Measurement Total 

Coordinates Timing Modeling 
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TABLE 1. RELATIVISTIC EFFECTS ON STATION POSITIONS AND GM^ 

ESTIMATED USING 6 TWO WEEK ARCS OF LAGEOS, SPACED TWO MONTHS APART IN 1980 
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TABLE 2. RELATIVISTIC EFFECTS ON YARLAS COORDINATES 

FOR ARC 5 


on nil coordinates in excess of half a cantiaeter, with a 
high degree of cancellation. The cancellation is aost 
apparent for height, for which there are four different 
contributions with aagnltudes in excess of 10 ca. The net 
effect on height is still only 1.6 ca. 

Table 3 suaaarises the individual arc siaulatlons for 
station coordinate estiaation, including the effect of using 
different values of GM^ for the relativistic and Newtonian 
data reductions. The last coluan in this table shows the 
arithaetic average of the effects on dX, 66, and 6H for each 
station. This coluan can be compared with the last column 
of Table 1, with which there should be good overall 
agreement. (The agreement should not be exact due to the 
slightly different weighting of each arc in the common 
estimation of station coordinates.) The rms agreement is 
.09 cm and the maximum difference is 0.2 cm, so the average 
station position change Is quite consistent with the common 
parameter solution. 

Examination of the month to month variations in Table 
3, however, shows that they can be at the several centimeter 
level for all three coordinates. The largest month to month 
change is in the YARLAS longitude, which shows a 6.6 cm 
variation between January 1980 and March 1980. Although not 
enough arcs have been simulated to obtain a clear picture of 
the variation during the year, it is expected that the 
estimated station coordinates would show smooth month to 
month (or arc to arc) variations. Due to the motion of the 
LAGEOS node, the effects would not be expected to have the 
same pattern the following year. 
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• 1 o Arc 3 Arc 4 Arc 5 Arc 6 

Station (o5n^980) (MarS 1980) (May 1980) (July 1980) (Sept 1980) (Mov 1980) 
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One other eepect of the elnuletlone auet elso be 
coneideredt the extent to which the relativletic effects 
would be expected to produce eyeteaetlc effects In 
aessureaent resldusls. Figures 1 snd 2 show the effects of 
the two aost slgnificsnt perturbstlons on aessureaent 
resldusls for two particular passes selected from the 6 arcs 
with individual arc station sdjustaents. Figure 1 shows 
residual effects for an ARELAS pass in Noveaber 1980, with 
the force aodel snd aessureaent effects tending to cancel 
throughout aost of the pass. There is still a net effect of 
over 7 ca during part of the pass* In Figure 2, which shows 
residual effects for a WETZEL pass in September 1980, the 
situation Is much worse. During aost of the pass, the 
measurement effects and force model effects have the same 
sign, and their sum has values as large as -20 ca. 

The Interpretation of the total relativistic curves 
in Figures 1 and 2 is as follows. If the data set for a two 
week period is first processed using the Newtonian theory, 
one set of residuals will be obtained. Applying the set of 
relativistic perturbations will then produce changes in the 
residuals as shown by the curves with solid dots in Figures 
1 and 2. If there were no other error sources, and the 
relativistic perturbations applied are correct (to within 
centimeter or so effects), then the residuals in the 
relativistic data reductions should be effectively sero. 
Thus, there must have been modeling errors (due to not 
including relativistic effects) in the Newtonian data 
reductions which produced residuals up to 20 cm or more. 
Such amplitudes far exceed the magnitudes of the station 
changes in Table 3. However, Table 3 is based on an average 






OmCMNALPAOEM 
Of POOR QOAUTY 



• fftet ov«r a l«rg« nuabar of paaaaa ovar a 2 waak parlod 
and 20 CB la a Baxlnuti raaidual affact on ona paaa, ao tha 
two raaulta ara not nacaaaarily ineonalatant. Mavarthalaaa , 
a 20 ca raaidual la a aufficlantly larga affact that raa 
diffaraneaa batwaan ralatlvlstlc data raductiona and non- 
ralativlatlc data reductlona ahould ohow algnlflcant 
diffaraneaa* 

It ahould be eaphaaicd that tha aiaulatlona dlacuaaad 
above have conaidarad only two weak arc langtha and have 
bean baaed on siaulatad data and not data actually taken 
during 1980. Although it la believed that the results are 
qualitatively valid for other arc lengths and for the use of 
actual data, this has not been verified. 
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SECTION 8.0 

SUMMARY AND CONCLUSIONS 


Algorlthaa for th« procoaaing of orbital tracking data 
baaad on tha ganaral theory of relativity hava baan 
praaantad. aaaunlng tha axlatanca of a conputar progran 
(auch aa Gaodyn) ImpIanantlAg Nawtonlan thaory. Laaar 
ranging data, auch aa that currently being taken tracking 
tha LAGEOS aatalllta. waa tha baalc ■aaauraaant type 
conaldarad • 

Tha dlffarancaa batwaan tha ralatlvlatlc thaory and 
tha Nawtonlan thaory can be briefly suanarlcad aa followa: 

1. Tha coordinate ayataa uaad ahould. at laaat 
lapllcltly. be cantered at (and aovlng with) tha 
aolar ayataa barycantar. alnca thla la the 
coordinate ayataa uaad In deriving tha 
ralatlvlatlc equatlona of aotlon of the aatalllta 
and planetary bodlea. 

2. Tha dlffarancaa between tha atoalc r.laa, kept at 
tha tracking aitea. and coordinate ti\,a. uaad In 
tha equatlona of aotlon. ahould be accounted for. 

3. The ralatlvlatlc equatlona of notion auat be 
laplaaantad . 

4. Tha light tine aolutlon appropriate to tha nobody 
aetrlc auat be uaad In relating the naaaured 
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round trip tiaes to the coordinetee of the 
eatelllte end tracking atation, 

5. Station coordlnatea uaad In the orbital solution 
must taka account of the velocity of the stations 
relative to the Integration coordinate systea. 

Of theae areas, accounting for atoalc tlaa-coordlnate time 
differences Is the closest to being negligible for LAGEOS 
data reductions, but It should be Included for 
coapleteness Item 5, the station coordinate transf oraatlon , 
Is probably In greatest need of theoretical analysis, since 
the transformation proposed has not been rigorously derived. 


Simulations show that use of the relativistic theory 
produces station position estimations which can vary at 
least several centimeters during the course of a year from 
estimates based on Newtonian theory. Averaging over a year, 
the differences are considerably reduced. The simulations 
also show differences In data reduction residuals (between 
the relativistic and Newtonian reductions) up to 20 cm In 
two week arcs in which all station positions are estimated. 


It should be noted that orbital data reduction 
programs developed for processing tracking data for earth 
orbiting satellites have traditionally been developed In 
earth-centered coordinate systems, and no change In this 
procedure Is recommended. However, the coordinate system 
used In the theoretical derivation of the relativistic 
equations of motion and light time solution is a solar 
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systen barycent«r coordinat* ayataa* Tha Intagration 
procaaa than producas tha coordinataa of tha aatalllta In 
tha barycantar coordlnata aystaa, ainua tha coordinataa of 
tha earth In tha barycantar coordlnata ayataa. Thia 
procadura allowa Incraaaad praciaion In aatalllta poaltlon 
and velocity coordinataa relative to tha earth. But aatel- 
lite velocitlea of 6-7 ka/aac relative to tha earth ahould 
not be confuaed with the velocity of the aatalllta - and 
earth - In the barycantar coordlnata ayataa^ which la of tha 
order of 30 ka/aec. For such a velocity, (V/c) « 10 , and 
this scale factor applied to characteristic distances for 
LAGEOS ( ~ 6 X 10^ a) corresponds to 6 ca. Relativistic 
affects of at least several centlaeters are thus to be 
expected on the basis of the velocities Involved, as well as 
froa gravitational effects. 

One modification (Section 6) has been developed for 
range measurement calculations to correct for an 
approximation originally made with the orbital program 
operating In an earth-centered coordinate system moving with 
the earth. It Is, of course, possible that other 
approximations are also no longer valid for the barycenter 
coordinate system. 

The ultimate test of the validity of a theory and Its 
Implementation Is for It to produce better data fits and 
more consistent estimates of physical parameters. This Is 
the obvious next step In the application of the Einstein 
gravitation theory to the reduction of laser ranging data to 
LAGEOS. 
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